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Abstract. The separability of the spatial modes of a charged particle in a Penning 
trap in the presence of an environment is studied by means of the positive partial 
transpose (PPT) criterion. Assuming a weak Markovian environment, described by 
linear Lindblad operators, our results strongly suggest that the environmental coupling 
of the axial and cyclotron degrees of freedom does not lead to entanglement at 
experimentally realistic temperatures. We therefore argue that, apart from unavoidable 
decoherence, the presence of such an environment does not alter the effectiveness of 
recently suggested quantum information protocols in Penning traps, which are based 
on the combination of a spatial mode with the spin of the particle. 



Submitted to: J. Phys. B: At. Mol. Phys. 

PACS numbers: 03.65.Ud, 03.67.Mn, 03.65.Yz, 37.10.Ty 

1. Introduction 

The possibility to confine a single charged particle in a Penning trap [TJ was 
initially widely exploited for high precision measurements of ion masses j2j[3] and the 
determination of fundamental constants [1H6]. More recently, however, Penning traps 
also began to attract the attention of the quantum information community, as a potential 
candidate for quantum computer building blocks [THH] . Indeed, the external control 
and cooling possibilities combined with the achievable high precision certainly present 
considerable advantages. It was shown earlier that quantum information protocols in a 
Penning trap are, at least theoretically, possible; The schemes known so far are based 
on information storage in either the axial or the cyclotron degree of freedom and the 
spin of the particle [TQrlT2] . Hence, as in any potential qubit implementation, one has 
to address the effects of an environment on the device. While decoherence is to some 
extent unavoidable, in the particular case of a Penning trap the environment may have 
yet another undesirable effect, as it allows for coupling of degrees of freedom which 
would not occur otherwise. More specific, the axial and the cyclotron motion are no 
longer uncoupled if an environment is present, and consequently the separability of the 
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latter cannot be taken for granted. It is the purpose of the present work to examine 
whether environment-assisted entanglement of these modes is likely to emerge. 

The intended analysis essentially requires a model which incorporates environmen- 
tal effects into the dynamics and a measure for separability (or entanglement) which 
can be extracted from such a model. For simplicity, we will assume a Markovian envi- 
ronment, and treat it by means of a standard master equation of Lindblad type [13J. At 
this point it should be mentioned that generally both a Markovian and a non-Markovian 
bath can (but not necessarily do) lead to environment-assisted entanglement [HHlH], 
so that this particular choice should not rule out an outcome a priori. Assuming linear 
coupling to the environment, the equations of motion for a charged particle in a Pen- 
ning remain analytically solvable. Furthermore, since the Penning trap Hamiltonian is 
at most quadratic in the canonical variables, initially Gaussian states remain Gaussian 
for all times. For Penning traps, such class of Gaussian coherent states was recently 
derived [Hl[20]. This makes it possible to establish a continuous variable separability 
criterion for the modes which is completely determined by the covariance matrix, as the 
general PPT criterion greatly simplifies for a system of 1 vs. N — 1 symmetric Gaussian 
modes. We will first describe the environment model in Section |2} from which the time 
evolution of the covariance matrix is obtained. We will then introduce a separability 
criterion which is extractable from the covariance matrix in Section |3j The results and 
conclusions will be presented in Section HI The work is summarized in Section |5j 



2. Phenomenological modelling of the environment 

The characteristic frequencies of a charged particle with mass m and charge q in a 
Penning trap are determined by the trap parameters (magnetic field B, trap dimension 
d and electrode voltage U) as follows: 

bj c = - — (cycltotron), uj z = \ (axial), (1) 

m V ma 2 

and the radial frequency is defined as 

The corresponding Hamiltonian can be written as 

H = H ± + H ZJ (3) 



with 



H± = 2m + -j(xp y - yPx) + ~y~ \ x + y ) ' ( 4 ) 
Hz =pL + I mw v, (5) 

2m 2 z ' V ; 

where the additional spin term is omitted, since the spin motion is completely separable 

from the dynamics and does not affect the calculations. We introduce the phase space 

vector 

R= (x,p x ,y,p y} z,p z ) T (6) 



Penning trap robustness against entanglement 



3 



and the covariance matrix of the system 



(7^ = -Tr \p(RiRj + RjRi)] - T^pR^pR,), i, j = 1, . . . , 6 (7) 

where p is the reduced density operator of the system. Its time evolution in the presence 
of a weak Markovian environment can be described by a master equation of the form 
dp i 

dt ~ ~h v " iri ' 2h 

J 

with Lindblad operators Vj linear in the operators pk and Xk 



(8) 



V j = Y,( a *jkVk + b* k x k ) 



J 



,6 



(9) 
and 



k=i it=i 
where a,-&, bj%. are complex coefficients. Models of this type are well known 
such an approach was adopted earlier to study damping and decoherence behavior in 
Penning traps [24J. However, unlike therein, here we will not simplify the picture by 
demanding that the damping in a certain degree of freedom may not affect the dynamics 
of another degree of freedom. By transforming the master equation to the Heisenberg 
picture, one obtains the following time evolution of the covariance matrix: 

a(t) = exp(At)(<r(0) - r)(exp(At)) T + T, (10) 

where cr(0) is the covariance matrix of the initially prepared state. The 6x6 matrices A 
and T are completely determined by the Lindblad coefficients a^, bj^ and the parameters 
of the Hamilton operator. As suggested in [22], it is convenient to introduce the vectors 



with the scalar product 



(b 



k = 1,2,3, 



<f,g> = 

i=l 

and the phenomenological constants (k, I = 1, 2, 3) 
an = -aik = -Im ((a k , ai)) , 
Pki = -Pik = -Im ((b k , bi)) , 
X M = -Im((a k ,bi)) , 



D 



D 



X[X k 



-Re((a k ,ai)) 



D 



pkPi 



D 



pipi,- 



a, 



'x k pi — ^pix k — 

The matrix A is then given by 



A 



Re((b k ,bi», 
-He^bi)). 



(11) 
(12) 

(13a) 
(136) 
(13c) 

(i3d) 

(13e) 
(13/) 



/ 


-An 


1/m 


-A12 - uj c /2 


-"12 


— A13 


-"13 






— muj\ 


-An 


P12 


-A21 - w c /2 


/?13 


-A31 






-A 2 i + co c /2 


«12 


— A22 


1/m 


— A23 


-"23 






-Pw 


-A12 + u c /2 


— mbj\ 


— A22 


/?23 


— A32 






— A31 


"13 


— A32 


"23 


— A33 


1/m 




V 


~Pl3 


— A13 


-fa 


— A23 


—moj 2 z 


— A33 


) 



(14) 
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and the symmetrical matrix T is the solution of the linear equation system [22] 

Ar + TA T = -2D. (15) 

Here, D denotes the symmetric 6x6 diffusion matrix, which is defined as Dij = Dr.^. 
It is also important to note that the phenomenological constants defined above cannot 
be chosen arbitrarily but have to satisfy certain conditions in order to preserve the 
positivity of the density matrix and the uncertainty relation. These restrictions are 
summarized in the appendix. 



3. Separability criterion for Gaussian states 

A necessary condition for the separability of a bipartite system was formulated in the 
1990s by Peres [25] and Horodecki et al |26j . They have proven that this information is 
directly related to the positivity of the partial transpose (PPT) of the density matrix 
with respect to one of the partitions. Although this is a powerful criterion, it is not 
always easliy applicable in practice in its original form. However, Gaussian multimode 
states have the unique feature that the complete information about the state can 
be extracted from the corresponding first two moments (i.e. the expectation values 
and the covariance matrix), without the necessity to deal explicitly with the density 
matrix. For practical purposes, this property is very valuable; In particular, it can 
be shown that the full uncertainty relation as well as conditions for separability or 
entanglement are extractable from the covariance matrix alone. Recent reviews with 
focus on entanglement in continuous variable systems can be found e.g. in [27,28j. 
Below, we only briefly summarize the most crucial aspects relevant for the present 
study. 

(i) Uncertainty relation 

Let us consider a system of A" Gaussian modes, with the phase space vector 
H N = (xi,pi, ..,xn,Pn)- We define the symplectic matrix Q as 

(16) 

Then the covariance matrix a, defined as in (JTj), has to satisfy the following so 
called bona fide condition in order to preserve the full uncertainty relation [29J: 

a + ^Q>0. (17) 

In other words, only covariance matrices satisfying the above condition are 
physically meaningful. For the particular case of a Penning trap, some related 
aspects were discussed in [30] . 

(ii) Separability 

Given a covariance matrix a that satisfies the bona fide condition, the following 
necessary and sufficient condition for separability holds for a system of m + n 
bisymmetric modes [3TH33] : 

a + ^tt>0. (18) 



n 



N 



Rf,R k 
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Here, a is the covariance matrix of the partially transposed state with respect to 
either of the partitions. The partial transposition in phase space corresponds to an 
operation which switches the sign of the momenta belonging to a partition, hence 
a is given by 

a = WaW, W = ®Ts z © hn, (19) 

where s z denotes the Pauli ^-matrix and i2 n is the 2n x In unity matrix. We 
emphasize that this holds, in particular, for the special case m — N — 1, n = 1 
which will be considered below. 

Given the above criteria, the time evolution of the covariance matrix ([TO]) thus provides 
us with the information whether environment-assisted entanglement occurs at any given 
time. This is the case whenever the condition (fl8|) is violated. 



4. Calculations, results and conclusions 



4-1. Preliminary remarks 



Following the notation from Eqs. ffT6l) - (fT9|) . the Penning trap is a particular case with 
N = 3, m = 2 and n = 1, since the cyclotron motion is completely symmetric in the x 
and y modes. The partial transposition with respect to the axial mode therefore becomes 
W = diag(l, —1, 1, —1, 1, 1) and is applied to the covariance matrix defined in ©. Since 
the matrix a + |f2 is Hermitean, it is sufficient to compute its lowest eigenvalue e in 
order to monitor the (violation of) non-negativity. The initial covariance matrix (7 > 
which enters the time evolution ( ITUl) . was calculated for the Penning trap coherent 
states in [T9|l34"]. The only remaining problem is to find an appropriate choice for 
the phenomenological constants (I13aj) -( T3^ ) , which is indeed one of the major obstacles 



when working with master equations of the form (jSJ). In the present case, the number 
of these parameters can be reduced for symmetry reasons; Furthermore, to some of the 
remaining constants it is rather straightforward to assign a physical meaning. Still, 
the phenomenology is not fully resolved. In what follows, we discuss this problem and 
suggest a Monte-Carlo-like approach to circumvent it. 

(i) Symmetry considerations and damping rates 

To illustrate the meaning of the elements of the matrix A defined in (|14p . we note 
that the equations of motion for the first moments (i.e. expectation values) derived 
from the master equation read 

A(R)( t ) = A (R)(t), (20) 
with the solution 

(R)(t) = exp(At)(R)(0). (21) 

We thus see that for k ^ I the parameters A/^ couple the coordinates Xk and xi and 
the momenta pk and pi, while the parameters a^i, (3ki couple the coordinate Xk with 
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the momentum p\. Since the system at hand is completely symmetric in x and y, 
we can reduce the number of parameters by setting 

A12 = A2I, A13 = A23, A31 = A32, 

«13 = «23, 013 = 023- (22) 

Moreover, for k = I the constants Xkk are identified as damping rates for the mode 
k, and, again for symmetry reasons, we can set An = A 2 2- 

(ii) Diffusion coefficients and temperature 

For the same symmetry reasons as above, the number of independent diffusion 
coefficients fll3dl) -( 13j) can be reduced by setting 

D X X Dyyi D xz Dy Z1 Dp x p x Dpyp^y 

Dp x p z — D PyPz , D xpy — D yPx , D xpx — D yPy , (23) 
D xpz — Dyp zl D zpx — D zpy . 

In addition, the temperature T of the environment can be incorporated by the 
following choice for the diagonal elements of the diffusion matrix: 

D xx = D yy = coth ( , D ZZ = ^- coth ( , (24) 



h\ 11 mu± I Hu± \ h\ 33 mu z ( hu z 
D„ „ = D v „ = coth , D v v = coth 



2k B T ) 1 

k B being the Boltzmann constant. This can be seen as a multidimensional extension 
of diffusion coefficients corresponding to an asymptotic Gibbs state. The choice is 
well-known in one-dimensional models of such type [2T | [23 |[35ll36] . 
(iii) Monte-Carlo approach 

Since the exact form of the environment is not known a priori, the computation 
of the lowest eigenvalue e of a + iQ/2 as a function of time for some particular 
choice of the remaining independent phenomenological constants would not really 
answer the question whether environment-assisted entanglement emerges; Instead, 
a scheme is implemented where a large number of trajectories e(t) is computed 
with the independent environmental constants being chosen randomly for each 
trajectory. The random generation is, however, restricted by the Born-Markov 
condition underlying the master equation (jSJ). This is taken into account when 
choosing the intervals from which the random variables are selected. In addition, 
the constraints listed in the appendix have to be satisfied for each trajectory. If 
this is not the case, the trajectory is discarded. The same applies if a violation of 
the bona fide condition occurs. The remaining physically meaningful trajectories 
e(i) obtained in this way allow us an insight into the entanglement dynamics and 
its dependence on the bath temperature. 
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4-2. Exemplary calculations and discussion 

As an example, we consider a proton in a Penning trap with the following parameters: 



The values are taken from Table II in pp. From here on, a system of units is used 
in which we set h, the proton mass and u z equal to unity. The positive damping rates 
^11 = ^22 and A33 are generated as uncorrelated, uniformly distributed random numbers 
in the interval [10 -2 ,10 -1 ], thus obeying the Markovian assumption Xkk w c , oo z 
while the coupling constants \m, ocm, (3m with k ^ I were assumed to be even weaker 
than the direct coupling terms and generated in the same manner from the manifold 
[1(T 3 , 1(T 2 ] U [-10~ 2 , -1(T 3 ]. To fulfill the conditions ( TO) we have for the off-diagonal 
elements of the diffusion matrix: 



where £_R fc R; £ [ — 1> 1] are a se t of random numbers. 

Figure [TJ shows the time evolution of the lowest eigenvalue e of a + iQ/2, computed 
for three different temperatures (10 mK, 0.1 K and 1 K). Each of the plots contains 
1000 trajectories (less the discarded non-physical ones). The diagrams were dissected 
into single bins, where the coloring indicates the number of trajectories passing 
through a given bin. We note that qualitatively the results are nearly identical, 
while quantitatively e scales linearly with temperature. The linear scaling is a typical 
sign of the high-temperature limit, as one obtains from the Taylor expansion of the 
diffusion coefficients ( 1241 coth. {Tiuj / {2k bT)) ~ T for ^> Tiuj. For the temperatures 

considered, the latter condition does indeed hold, and hence one would expect the 
detrimental effect of thermalization on entanglement to be dominant compared to the 
weak environmental coupling of the modes. The results shown in Figure [1] confirm the 
expectation - from the diagrams it is evident that, even if possible, in the studied 
case environment-assisted entanglement is very unlikely, as for all the trajectories 
there is not a single event of e becoming negative at any time. We note by passing 
that the same observation holds also if the number of trajectories is increased by 
3-4 orders of magnitude, which just leads to a rescaling of the plots, leaving the 
number of entanglement events unchanged equal to zero. Thus, at the typical operation 
temperatures, the environment is much more likely to destroy rather than create 
entanglement in a Penning trap. In order to corroborate this hypothesis, we also 
considered the very low temperature case T = 1 mK, which, to the best of our knowledge, 
is not yet experimentally feasible for a Penning trap. In this case the high-temperature 
limit breaks down and consequently the thermalization is less severe, thus possibly 



u c » 484 MHz, u z « 63.2 MHz. 



(25) 




(26) 
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T=10 rnK T=0.1 K T=1 K 




Figure 1. (Color online) Monte-Carlo trajectories of the lowest eigenvalue e of the 
matrix a + iil/2, where a is the partially transposed covariance matrix, shown as a 
function of time for three different temperatures. The coloring indicates the number 
of trajectories passing through a given bin; We particularly emphasize the absence of 
negative values of e (all trajectories start at e(t = 0) = 0), which persists even for 10 6 
computed trajectories (not shown here). The time unit in the chosen unit system is 
ps 15.8 ns. 

allowing for environment-assisted entanglement. Indeed, further numerical simulations 
support this conclusion. Figure [2] shows 10 7 trajectories (less the discarded non-physical 
ones), computed in the same manner as those displayed in Figured] but at T = 1 mK. 
One clearly observes that the non-negativity of a+ifl/2 is violated for some trajectories, 
which is a signature of entanglement. To ensure that the entanglement is solely due to 
the temperature decrease and not a particular constellation of random numbers, we 
also checked that the sets of random numbers that give rise to entanglement at low 
temperatures do not do so at higher temperatures. Therefore, we can conclude that 
our results strenghten the assumption of separability of the cyclotron and axial mode 
in quantum infromation related applications of Penning traps, at least in the case of a 
weak Markovian environment at realistic temperatures. 

5. Summary 

We investigated the possibility of environment-assisted entanglement of the spatial 
modes of a single charged particle in a Penning trap. The system-environment 
interaction was modelled by a Markovian master equation, with Lindblad operators 
preserving the Gaussian form of an initially prepared Penning trap coherent state. This 
allowed us to apply a specific practical form of the PPT separability criterion extractable 
from the covariance matrix, the time evolution of which is analytically solvable in the 
framework of the adopted model. The generally unknown phenomenological parameters 
of the master equation were treated by means of a Monte-Carlo-like simulation. We 




Figure 2. (Color online) Same plot as in Figure 1 but at T = 1 mK, shown for 
10 7 trajectories on a logarithmic scale (however, for a better visualisation the shading 
(color) of the bins containing zero trajectories is set to zero instead of — oo). At this 
low temperature, the emergence of entanglement is clearly visible. 



found that the environment is very unlikely to cause entanglement in the considered 
case, implying robust separability of spatial Penning trap modes. By comparing the 
time evolution for different temperatures, this behavior should most likely be attributed 
to the generally destructive effect of thermalization on entanglement. This is supported 
by entanglement signatures occuring at unrealisticly low temperatures. Thermalization 
thus turns out to be much more significant than the weak coupling of the modes caused 
by environmental scattering. 

Note added: After finishing the manuscript, we became aware of two very recent 
papers which should be mentioned in connection with the present work. The first paper 
by Hamdouni [37] contains an explicit analytical derivation of multidimensional diffusion 
coefficients for the same environment model as used here. This result links the diffusion 
coefficients to the elements of A (Eq. ( fT4]) ). which in principle allows us to restrict the 
Monte-Carlo approach solely to the phenomenological constants given in ( 1221) . The 
second paper by Anastopoulos et al [38] addresses the generalized uncertainty relations 
and generation of entanglement in quantum Browninan motion using Wigner propagator 
techniques; our results qualitatively agree with several findings of their work. 
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Appendix A. Restrictions on the phenomenological constants 

The constraints are a straightforward three-dimensional extension of the relations known 
for the two-dimensional case [22J. The parameters have to be chosen such that the 
principal minors of the following matrix are non-negative: 



/ 




n iha 12 
J^xy 2 


^xz 2 


n ifiAn 
-^xpic 2 


n iRAia 
-^Pj, 2 


D X p z 


2 


\ 




n 1 ihai2 
J^xy T 2 


Dyy 


n i^«23 
-^yz 2 


n iftAai 

•^HPs 2 


n «?iA22 
^S/Pb 2 


~Dyp z — 


iH\23 
2 






n _|_ *ft<*l3 

-LSXZ ~ 2 


1 2 




n i?iA3i 

^zp x 2 


n i?iA32 
-^zPy 2 


~D zPz — 


iS,A33 
2 






n 1 

-^xpx T 2 


■^J/Pi "i" 2 




Dp x p x 


-^P^P:, 2 


Dp x p z ~ 


2 






n 1 i^i2 
^xpy T 2 


n 1 i^A 2 2 
^S/Pw ~ r 2 


n i i?iA 32 

-t^ZPy "I" 2 


n 1 ihpi2 

^PxPy ' 2 


A>3/P» 


DpyPz ~ 


ih/323 
2 




V 


n _|_ ihx ^ 


n 1 «?iA23 

up? T 2 


^zp z T 2 


J^PxPz ~ 2 


^PyPz ' 2 


Dp z p 


Z 


) 



In particular, the definitions fl 1 3 - ( fl3^ ) imply, because of the Cauchy-Schwarz 
inequality, the following constraints (k, I = 1,2,3): 

h 2 al 



D D - D > 

^ X k X k ^ XlXl ^x k x t 



HI 



h 2 \ 2 



D D — D z > - ' kl 

J^X k X k ^'p l p l ^XkPi — ^ • 

Finally, we note that the system only approaches an asymptotic state if the matrix A 
defined in (TH|) has no eigenvalues with positive real parts. In the illustrative Monte 
Carlo calculations, all trajectories e(i) computed with a set of parameters violating any 
of the above conditions were discarded. 
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